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Abstract 

The Buneman instability occurring when an electron population 
\ is drifting with respect to the ions is analyzed in the quantum lin- 

CN ' ear and nonlinear regimes. The one-dimensional low-frequency and 

collisional model of Shokri and Niknam [Phys. Plasmas, 12 (2005) 
062110] is revisited introducing the Bohm potential term in the mo 
^ ' mentum equation. The linear regime is investigated analytically, and 

^ ■ quantum effects result in a reduction of the instability. The nonlinear 

regime is then assessed both numerically and analytically, and pure 
quantum density oscillations are found to appear during the late evo- 
lution of the instability. 



The Buneman instability pQ is a basic instability process in classical plas- 
mas. It occurs in beam/plasma systems when there is a significant drift be- 
tween electrons and ions. It is frequently referred to as the "Farley-Buneman" 
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instability because it was almost simultaneously discovered by Farley PJ. In 
view of the very basic "set-up" it needs to be triggered, this instability has 
been found to play a role in many physical scenarios. In space physics and 
geophysics, this instability has been invoked in the earth ionosphere [2J or in 
the solar chromosphere [3|. In double-layer and collisionless shock physics, 
the same instability has been found responsible, under certain circumstances, 
of the very formation of this kind of structures [I]. In situations where an 
electron beam enters a plasma, like in the fast ignition scenario for inertial 
fusion [5J , the electronic return current prompted in has been found to gener- 
ate exponentially growing Buneman modes through its interaction with the 
plasma ions [6, 7]. 

On the other hand, presently quantum plasmas are attracting much at- 
tention from a variety of reasons, since they provide a typical state of ionized 
matter under large densities and/or low temperatures. For instance, quan- 
tum effects in plasmas are relevant in ultra-small semiconductor devices, 
metal clusters, intense laser-solid interaction experiments and compact as- 
trophysical objects like neutron stars and white dwarfs (see JHUH] for reviews). 
Moreover, X-ray Thomson scattering techniques in dense plasmas have been 
used [TUj to verify the signature of quantum diffraction effects in the disper- 
sion relation for electrostatic waves. Also, in the near future the development 
of coherent brilliant X-ray radiation sources [TTJ and keV free electron lasers 
[T2] will provide experimental access to the quantum nature of plasmas under 
extreme conditions. Some of the most recent developments in the field are 
the analysis of wave breaking in quantum plasmas [13], the characteristics 
of bounded quantum plasmas including electron exchange-correlation effects 
|14j . the development of a quantum single- wave theory for nonlinear coher- 
ent structures in quantum plasmas [15] . the discussion of waves in quantum 
dusty plasmas [T6], the prediction of a fundamental size limit for plasmonic 
devices due to the quantum broadening of the transition layer [17] , as well 
as the inclusion of spin [IS] and relativistic [TP] effects in quantum plasma 
modeling. Finally, we note the usefulness of quantum plasma techniques 
to other, closely related problems, such as the treatment of nonlinear wave 
propagation in gravitating Bose-Einstein condensates [20] , 

The aim of the present work is to discuss the quantum analog of the 
Buneman instability. Our approach is based on the quantum hydrodynamic 
model for plasmas, which has proven to be very useful for the understanding 
of nonlinear problems in quantum Coulomb systems [21] . More specifically, 
we consider linear and nonlinear low-frequency waves in a collisional electron- 
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ion plasma, extending the model by Shokri and Niknam [22] by means of the 
inclusion of a quantum term (the so-called Bohm potential) associated to 
the wave nature of the quantum particles. As verified in the continuation, 
the Bohm potential has a stabilizing influence on the low-frequency linear 
Buneman instability. In addition, it eventually produces nonlinear oscilla- 
tory structures at the late stages of the instability, in sharp contrast to the 
monotonic character of the classical stationary states. 

We start with the two-species cold quantum hydrodynamic model [23J, 
taking into account the effect of collisions using simple relaxation terms, 
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InEqs. ([I])-© the electron (ion) number densities, v e> i the electron 

(ion) fluid velocities, E the electrostatic field, m (M) the electron (ion) mass, 
— e the electron charge, Ti the Planck constant divided by 2ir and £o the vac- 
uum permittivity. Moreover, v e ^ represent electron (ion) collision frequencies 
with neutrals. For simplicity in this work only one spatial dimension is con- 
sidered. Quantum effects are included in the force equation for electrons by 
means of the ~ h 2 term, the so-called Bohm potential. Due to m/M <^ 1, 
no quantum terms are needed in the ion force equation. 

Linearizing the model around the homogeneous equilibrium 

n e>i = n , v e = , Vi = — — , E = E , (6) 

mv e M Vi 

where E is an external DC electric field. Supposing perturbations ~ exp(i[kx- 
ut}) with wavenumber k and wave frequency u, the result in the reference 
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frame of the drifting ions is 

1 _ < 



(u — kvo)(u — kvo + We) — ft /c 4 /(4m 2 ) 

-0, (7) 



id [id + IVi) 

where uj pe = (n e 2 / (m£ )) 1//2 and u P i = (n e 2 / (Meo)) 1 ^ 2 are resp. the elec- 
tron and ion plasma frequencies and where 

= ~eE (— + -M (8) 
\mv e MPiJ 

is the relative electron-ion equilibrium drift velocity. 

In the very low frequency range u < i/j < kv , v e <C kvo, the dispersion 
relation reduces to 

= iu 2 pi k 2 (v 2 - h 2 k 2 /{Am 2 )) 
U v t id 2 e - k 2 v 2 Q + Ti 2 k A /{Am 2 ) ' 1 ' 

This mode is unstable (Im(w) > 0) provided 

<>W -^>*> ( 10 ) 

otherwise it is damped. Small wavelengths such that Ti 2 k 2 > 4m 2 t>Q are 
automatically stable, due to the quantum effects. 

Assuming a large ion- neutral collision frequency where z/, 3> u P i, we can 
suppose a slow temporal dynamics. Moreover, from Eq. (JTU1) we have vo/u) pe 
as a natural choice of spatial scale for the development of the instability, at 
least for not very large quantum effects. Therefore, we consider the following 
rescaling, 

toLt u pe x v e 

x-*^—, v e ^—, (11) 



pi 



Vi v v 

MviVi n ei eE 

Vi — > , n e i — >■ — , hi — > . 

mujpeVo n mujpeVo 

For simplicity using the same symbols for original and transformed vari- 
ables, the rescaled system reads 

uj 2 i dn e d(n e v, 



ViUpe dt dx 



+ \ 7 = 0, (12) 
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^ + 8^) = Q (13) 

u 2 pi dv e dv e v e 

— 7TT + V e~a~ = ~ E V z ( 14 ) 

ViOJpe Ul OX U! pe 

h 2 d (d 2 ^rjdx 2 ^ 

2 dx \ x /n, 
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dE 

dx 



where 



• r — } = E-Vi, (15) 
m-n e , (16) 



^ = 3 ( 17 ) 



•'ye. 
UIVq 

is a non-dimensional parameter measuring the relevance of the Bohm poten- 
tial. 

Provided 

Ui^>uJ P i, w pe >z/ e , (18) 



we obtain 



d{n e v e ) 
dx 



0, (19) 



(23) 



£ = u», (22) 

as 

where the remaining terms are assumed to be of the same order. The final 
equations are the same as Eqs. (2)-(6) of ref. [22], with the inclusion of the 
extra ~ H 2 contribution. The purpose of the present work is to investigate 
the role of this quantum term. 

We linearize Eqs. ( TE?j) - (l2"3"j) around n e = rii = l,v e — 1, Vi = 0, E = 0. 
Note that after rescaling the equilibrium ion velocity and electric field are 
higher-order terms, due to v e juj p( , <C 1. We get the dispersion relation 

k 2 (l-H 2 k 2 /A) , , 
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Figure 1: Growth rate 7 in Eq. ( 12"4"|) as a function of the wavenumber k, 
in the classical limit (H = 0). Dimensionless variables are used. Note the 
asymptote at k = 1. In addition, 7— >— lasA;— >oo. Instability is found for 
< k < 1. 

which is the same as Eq. fl9]), in terms of non-dimensional variables. Moreover 
in Eq. ( !2~4"1) we define 7, the imaginary part of the frequency. 

It is interesting to analyze the behavior of 7 according to the quantum 
parameter H. In the classical H = case, one has 7 = k 2 / (1 — k 2 ) and linear 
instability for k < 1. For the sake of comparison with the non- vanishing H 
case, in Fig. [T]we show the corresponding form of the classical linear insta- 
bility. The asymptote at k — 1 points for an explosive instability. However, 
this singularity is eventually regularized by nonlinear effects, as discussed in 

In the semiclassical < H < 1 case, the growth rate from Eq. (I24p has 
two asymptotes at kA,B defined by 

H 2 } 1 ' 2 ) , 

H 2 ] 1 ' 2 ) . (25) 

Moreover, one has 7 > for < k < kA or ks < k < kc, where 

k 2 c = k 2 A + k 2 B = ^- 2 . (26) 

By coincidence, one has formally the same instability condition as for the 
quantum two-stream instability described by a quantum Dawson model, see 
Eqs. (37)-(42) of ref. [23]. We note the instability of small wavelengths 
where ks < k < kc has no classical counterpart. The behavior of 7 as a 
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Figure 2: Growth rate 7 in Eq. f [2"4"l) as a function of the wavenumber k, in 
the semi-classical case (0 < H < 1). The value H = 0.5 and dimensionless 
variables were used. Note the asymptotes at k = fc^s. In addition, 7 — > — 1 
as fc — > 00. Instability is found for < k < and also in the small 
wavelength region k% < k < kc- 



function of the wavenumber in the semiclassical situation when < H < 1 
is shown in Fig. [2J Again, 7 — > — 1 as k — > 00. 

The case H = 1 is particular because then k^ = kg = a/2, so that the 
mid stable branch in Fig. [2] disappears. One still has explosive instability, at 
k = \/2. The perturbation is linearly stable for k > kc = 2. This is shown 
in Fig. 

When H > 1, the denominator in Eq. (1241) can be shown to be always 
positive, so that singularities are ruled out. Instability is found for k < 
kc = 2/if, with the most unstable wavenumber being k = y/2/H. The 
corresponding maximal growth rate is 7 max = ^/(H 2 — 1). We see that 
both the unstable k— region and the maximal growth rate shrinks to zero as 
H increases, which is a signature of the ultimate stabilizing nature of the 
quantum effects. The corresponding function 7 (A;) is shown in Fig. |H 

The overall situation can be visualized in Fig. 5, where the unstable 
region in (k 2 , H 2 ) space is shown. This is formally the same as Fig. 1 of 
ref. [23] on the quantum two-stream instability. On the other hand, for 
the low-frequency collision-dominated Buneman instability one has distinct 
behaviors of 7 (A;), according to the parameter H. Namely, one has the four 
classes shown in Figs. [l]to|H m contrast, the quantum two-stream instability 
exhibits no singularities of the growth rate at specific wavenumbers. 

The linear instabilities just described are eventually killed by nonlinear 
effects. In this context, Eqs. (fT§|) - (l2"5|) provide a convenient framework for 
nonlinear studies of the collision-dominated low-frequency quantum Bune- 
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Figure 3: Growth rate 7 in Eq. (1241) as a function of the wavenumber k, in the 
particular case H = 1. Dimensionless variables are used. Note the asymptote 
at k = Ica = ks = V%- In addition, 7—^—1 as k — >■ 00. Instability is found 
for < k < 2. 
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Figure 4: Growth rate 7 in Eq. (|24p as a function of the wavenumber k, in 
the strongly quantum (H > 1) case. The value H = 2 and dimensionless 
variables were used. Again, 7— >■ — 1 as -> 00. Instability is found for 
< k < k c = 2/H. 
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Figure 5: Stability diagram for the low-frequency collision-dominated quan- 
tum Buneman instability. The filled area is unstable. Lower and middle 
curves: resp. k\ and k 2 B as defined in Eq. ( 1231) . Upper curve: as defined 
in Eq. 



man instability. We define the new variables 

1 



V 



p = rii-n e , 



(27) 



giving resp. the electron fluid velocity and net charge density. Equations 
(TP9"j) - (l2"5j) can then be shown to reduce to 



1 d 2 
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~p d ( 
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2dx\ 



V 2 -H 



2 d 2 ^l/V/dx 2 



l/V 



Assuming quasineutrality (p = 0) and V — 1 at early times, we find the 
first three terms in Eq. f[2l?|) to be initially the more relevant. Linearizing 
them we get 

dV d 2 V H 2 d 4 V ^ 



dt 



dx 2 



4 dx 4 
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a quantum-modified diffusion equation with a negative diffusion coefficient. 
Growing in time solutions are easily found, with the corresponding instability 
described by the dispersion relation 

, 2 f H 2 k 2 \ , . 

w = ik 2 ( 1 - —J . (31) 

This is the quasineutral version of Eq. (j24l) . However, the instability in 
time is accompanied by a periodic structure in space, a feature not present 
in the classical case. Indeed, oscillatory in x solutions to Eq. (I3"U1) can be 
also readily constructed e.g. by separation of variables. In addition, again 
quantum effects provide the stabilization of the small wavelengths such that 
k> k c = 2/H. 

After the initial increase of the perturbation and enlargement of the den- 
sity gradient, the ~ p terms in Eq. (1291) become essential. To examine the 
stationary states of the model, we set all time-derivatives to be zero. A little 
algebra then shows that 



d 2 ( V 2 H 2 d 2 ^l/V/dx 2 \ 1_ 
dx 2 \ 2 2 JTJV I ~ V 



assuming symmetric solutions so that V"(0) = and excluding the case of 
identically vanishing electric fields. 

We start solving Eq. (1321) in the classical (H = 0) limit. Defining 

V 2 

K = - (33) 



one obtain the Newton-like equation 

d 2 K 



dx 2 J2K 



(34) 



Assuming n e (0) = 1 and a symmetric density profile so that n' e (0) = 0, one 
can integrate Eq. (1341) twice with K(0) = 1/2,K'(0) = 0. In terms of the 
electron fluid velocity V the result is 

(y-l)^(V + 2)=X, X = ljL, (35) 
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which is equivalent to a cubic equation for V with only one real root, namely 



V = -l + Q + l/Q, (36) 

where 

Q= ( 1±MH^ '\ (37 ) 

This is the classical nonlinear stationary Buneman solution, in full agreement 
with Ref. [22]. We note that the ion profile follows from V{ = E, with E 
given by Eq. f[2Tj) . as well as from m = p + l/V, with p given by Eq. f[28"j) . It 
turns out that rij = 0, which is a result from the low- frequency assumption. 
Indeed, the continuity equation for ions imply rtiE = cte., set to zero in view 
of the larger ion mass, which in turn imply a negligible stationary ion density. 
On the other hand, using Eq. (1361) one obtain an electric field E ~ — (6a;) 1 / 3 
for large \x\, arising from the electron fluid bunching. Also observe that 
more general solutions are possible to find, including ion density corrections 
as well as traveling wave forms. However, in this case the algebra becomes 
much more involved. 

Incidentally, we note that from Eq. fl3"6"|) we have 

n e (0) = l, <(0) = 0, <(0) = -l, <"(0) = 0, (38) 

a result to be used in the H ^ case. 

In the quantum case, we need to return to Eq. (I3"2|) . Defining 



A = ^T e = ^/T/V, (39) 
one obtain the fourth-order ordinary differential equation 



2 d 2 (d 2 A/dx 2 \ d 2 ( 1 



"**{— —)-^{a<) + 2A =°< (40) 

describing the final, stationary nonlinear stage of the density modulations. 
For the sake of comparison, we assume 

A(0) = 1 , A'{0) = , A"(0) = -1/2 , A'"(0) = , (41) 



which correspond to the same boundary conditions (138j) of the classical so- 
lution. In other other words, the classical solution is used to set also the 
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Figure 6: Stationary electron fluid velocity. Upper, left: classical exact so- 
lution from Eq. f)36p . The remaining comes from the numerical solution of 
Eq. (HOD. Upper, right: H 2 = 0.5. Bottom, left: H 2 = 1.0. Bottom, right: 
H 2 = 1.5. 

second and third-order derivatives of the density at x — 0, which are needed 
in the quantum case. 

Equation (HOI) can be numerically solved, yielding periodic modulations 
of the stationary velocity and density profiles. The amplitude of the modula- 
tions is seen to increase with the strength of the quantum effects, as apparent 
from Figs. [6] and [7J The emergence of new oscillatory structures is ubiqui- 
tous in quantum plasmas. In the present case, the ultimate role of the Bohm 
potential term in Eq. f[2~Tj) is a qualitative modification of basic equilibrium 
macroscopic properties like the electron fluid density and velocity. Similar 
oscillatory patterns of a pure quantum origin appear in the description of 
weak shocks in quantum plasmas [24J, of the quantum Harris sheet solution 
in magnetized quantum plasmas [25] and in undulations of the equilibrium 
Wigner function in quantum plasma weak turbulence [26J. 

Considering the extreme quantum case where only the ~ H 2 term is taken 
into account in Eq. (41) can be instructive for the physical interpretation 
of the influence of the Bohm potential. In this H 2 ^> 1 situation, one has 
A = cos(a;/v / 2) in view of the initial conditions in Eq. (42). This imply a 
density n e = A 2 oscillating with a wavelength A = a/2 tt, of the order of v /co pe 
using dimensional coordinates. The behavior so described is confirmed by 
numerical simulations, with corrections arising from the classical terms. For 
intermediate values of H, numerical analysis shows the wavelength displayed 
in Figs. El and [7J is actually not constant, but falls like ~ a; _3//4 instead. 

Turning attention to the quasineutral regime, we examine the nonlinear 
stationary states when p = 0. In this case, Eqs. ( l2"g|) - (T2"9|) reduce to Eq. 
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Figure 7: Stationary electron fluid density. Upper, left: classical exact so- 
lution from Eq. (I3T)|) . The remaining comes from the numerical solution of 
Eq. (Hfl). Upper, right: H 2 = 0.5. Bottom, left: H 2 = 1.0. Bottom, right: 
H 2 = 1.5. 



( 1281) only, the other one being redundant. After twice integrating using Eqs. 
(141jl . the model can be shown to be equivalent to the autonomous Pinney's 
[2"T] equation 

TT2 d 2 A / # 2 \ , 1 

*** + { 1 + T) A= #- < 42) 

Pinney's equation is endemic in nonlinear analysis and is well-known to be 
exactly solvable. This is specially true in the present autonomous case, where 
Eq. (l4"2"j) can be directly integrated twice. Assuming A(0) = 1,A'(0) = as 
before, the solution reads 



1 ( , H 2 „ r. zzrz: x 



A 2 = n e = ^r— ( 1 + — cos 2 

l + H 2 /2\ 2 



1 + H 2 /2 - 



(43) 



displaying quantum oscillations not existing in the classical case. Once again, 
the amplitude of the quantum oscillations in space increases with H, as shown 
in Figs. Owith the stationary electron fluid density. The bunching present in 
the non-quasineutral case is eliminated. Similar results apply to the electron 
fluid velocity. 

In conclusion, we have analyzed the low-frequency collisional quantum 
Buneman instability, both in the linear and nonlinear regimes through a 
one-dimensional model. Note that because this electrostatic instability is 
longitudinal, such a low dimensional analysis is relevant. The nonlinear evo- 
lution of the instability can be studied numerically and analytically, and 
results in pure quantum density oscillations. 
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Figure 8: Stationary electron fluid density in the quasineutral de- 
scribed by Eq. (H3l) . Upper, left: if 2 = 0. Upper, right: H 2 = 0.5. Bottom, 
left: H 2 = 1.0. Bottom, right: H 2 = 1.5. 
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